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We study the nonlinear Schrodinger equation with a periodic delta-function potential. This 
realizes a nonlinear Kronig-Penney model, with physical applications in the context of trapped 
Bose-Einstein condensate alkaly gases and in the transmission of signals in optical fibers. We find 
analytical solutions of zero-current Bloch states. Such wave-functions have the same periodicity of 
the potential, and, in the linear limit, reduce to the Bloch functions of the Kronig-Penney model. We 
also find new classes of solutions having a periodicity different from that of the external potential. We 
calculate the chemical potential of such states and compare it with the linear excitation spectrum. 
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Nonlinearity can deeply modify the Bloch theory of 
non-interacting atoms trapped in periodic potentials. 
Loop structures, energetic and dynamical instabilities, 
solitons and "generalized" Bloch states (i.e. states which 
do not share the same periodicity of the lattice), all 
arise in the context of a nonlinear Schrodinger (or Gross- 
Pitaevskii) equation. Applications span, for instance, the 
physics of dilute Bose-Einstein condensed gas trapped in 
optical lattices [1-9] or the propagation of signals in op- 
tical fibers [10] 

In this paper we study analitically the nonlinear 
Schrodinger equation with an external Kronig-Penney 
(KP) potential, given by a periodic array of delta- 
functions. The linear Schrodinger equation with the same 
potential has been solved quite early in the '30, playing 
a distinguished role as a model in metal's theory [11,12]. 

It is noticeable that also several properties of the non- 
linear Schrodinger equation, with the same KP external 
potential, can be derived analitically. The most interest- 
ing results, however, are related with the emergence of 
new properties which do not have a counterpart in the 
linear case. Stationary solutions of the GPE which do 
not reduce to any of the eigcnfunctions for a vanishing 
nonlinearity have been studied in [6] using a tight binding 
approximation, in and two- wells systems, in [13-17]. 

The mean- field model of a quasi- ID BEC trapped in 
a KP potential is governed by the following nonlinear 
Schrodinger (or Gross-Pitaevski) equation: 
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where p is the chemical potential and g the nonlinear 
coupling constant. The KP external potential is given by 
equispaced delta-functions: V(x) =p X)^L-oo 5{x—na), 
having a lattice constant a. Since the external potential 
has a step-like shape, it is useful to rewrite the GPE in 
hydrodynamic form. With if> (x) = \fp (x) exp [— i0 (x)] 
(and in dimensionless units), we have: 
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(/i, V), ap — » p and normalization J 1=c± dx p(x) = 1. Nq 
is the number of atoms in each well, and the integration 
constants a, [3 are fixed by the boundary conditions. In 
particolar, a has a simple physical meaning, being the 
current carried by the order parameter ij)(x): 



J = a 
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Bloch state. We derive a class of stationary solutions 
of Eq.s(2) having (i) the same periodicity of the external 
potential (Bloch states), and (ii) the linear limit: such 
states reduce to the well know solutions of the linear KP 
model when r\ — > 0. The Bloch theorem assures the possi- 
bility to write the complete set of eigenstates of the linear 
equation in the form: ip (x) = exp (iqx) f q (x), where q 
is the quasi-momentum of the state and f q (x) is peri- 
odic with the lattice constant, f q (x + a) = f q (x). It is 
immediate to extend the Bloch theorem and show that 
stationary states having the same periodicity of the po- 
tential exist also in the nonlinear case. However, there 
are new classes of states which do not have such proper- 
ties, see next section. 

Bloch nonlinear stationary states of the GPE with the 
KP external potential can be written in terms of the Ja- 
cobi Elliptical functions. With the same class of functions 
have been previously studied the exact stationary solu- 
tions of the GPE descibing a train of solitons [18], or the 
stationary solutions of the GPE in a double square-well 
potential [13-15]. To the best of our knowledge, however, 
the nonlinear Schrodinger equation with a periodic KP 
potential has not been studied so far. We notice that our 



results with the delta- function potential can be general- 
ized (although becomes rather cumbersome) to a regular 
array of square- wells, as will been shown elsewhere. 
Stationary solutions of the Eq.s(2) are: 



To solve Eq.(7) we need a further relation between the 
chemical potential p and K, provided by the normaliza- 
tion in (6). We obtain two conditions associated with the 
+1 and the —1 of Eq.(7), respectively: 
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and SN(u, n) is the Jacobian elliptic sine function. 
SN(u, n) has the desiderable property to give, in the 
linear limit, r\ = 0, p (x) = §^ cos 2 (Kx + 5) + 

8^- sin 2 (Kx + S), 9(» = axctaxi(8KaUn(Kx + S)/A), 
which are the exact solutions of the linear KP model. 

The three parameters A, K and 5 are fixed by impos- 
ing the continuity of the order parameter and the Bloch 
periodicity, as in the linear case. Because of the nonlin- 
earity, however, we need one more condition to fix the 
chemical potential: the normalization of the density. We 
therefore obtain four conditions: 
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-(l-n)K = 0(-) (9) 



where n = | (l — p/K 2 ^, tp 2 = K.(n) + ^ and ipi — 
fZ (n) — -y. £(u,n) is the incomplete elliptic integral of 
the second kind and K, is the complete elliptic integral of 
the second kind. Am (4>i, n) is the amplitude of 4>i- 

The coupled equations (7) and (9) can be solved graph- 
ically as shown in Fig.(l). The black lines are solutions of 
(7), while the color lines are solution of (9)) (the dashed 
and the continuum lines correspond to the +, — sign, re- 
spectively). The lines among the intersections between 
the black and the color lines (evidenciated by the color 
filled regions), give the values of (p, K) corresponding to 
the zero-current Bloch states. As expected, nonlinear- 
ity quite modifies the lower energy bands of the systems, 
while, for higher bands, differences are reduced. 
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In this paper we consider the special case of zero-current 
states, having a — (so that = const.). Using the 
elementary properties of the Jacobian Elliptical function 
[19] the first two equations of (6) give: 
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In the linear limit such states are at the top or at the 
bottom of the corresponding energy bands (which is not 
necessarly true in the nonlinear case when loop- like struc- 
tures appear in the excitation spectrum). Indeed, the 
Eq.(7) with n = (giving u = K 2 ) reduces to the well 
known relation 
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giving the band and the gap widths in the linear K-P 
model [12]. 
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FIG. 1. Graphical solution of the system of Eq.s(7) and 
Eq.s(9) for different values of nonlinearity. The allowed val- 
ues of the chemical potential p as a function of momentum 
K are given by the intersections between the black and the 
color lines, in the color filled regions 



Generalized Block states. In the previous section we 
have considered Bloch-like solutions, where the conden- 
sate wavefunction have the same periodicity of the po- 
tential. As already mentioned, the nonlinear interac- 
tion allows for stationary solutions which can have, in 
principle, any integer value period. In the following, 
we only consider solutions which are periodic every two 
sites: f q (x + 2) = f q (x). Notice that wave functions 
with q — 0,7r are symmetric, ip (x) = ip(x + 2), while 
are antisymmetric when q = n/2, ip (x) = —ip(x + 2). 
Therefore, f q= o(x) = f q=7r (x), and the respective chemi- 
cal potentials are equal. 

In the following, we consider as elementary cell two 
neighboring wells separated by the delta potential, with 
pi (x) and p r (x) being the densities in the "left" and 
"right" well, respectively. As in the previous section, we 
consider only states with zero current a = 0. The conti- 
nuity and the periodicity conditions give: 
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FIG. 2. Symmetric generalized Bloch states with no nodes 
(a), and with two nodes (b). 
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First, let's consider the solution which have nodes at the 
boundary of the elementary cell: 



Solutions with a node in x = 1 can be constructed from 
the previous conditions and imposing pi (x) = p r (x + 1), 
p r (x) = pi(x - 1). 



There is a different class of generalized Bloch states 
having nodes located outside the boundaries of the po- 
tential. The conditions are: 



Pl (0) = Pr (2) = 0. 
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Then we have 6i = (2Zi + 1) K (m), S 2 + 2K 2 = 
(2l 2 + 1) /C {n 2 )- Notice that the boundary conditions are 
equal to those of a double well, except for the condition 
on the first derivativies at the borders of the elemen- 
tary cell, which gives A{1 — n{) = B(l — n 2 ). Com- 
bining with the normalization, we arrive at A — B and 
K = Q. Notice that this relation excludes the existence 
of symmetry broken solutions which are instead obtained 
in a single double- well potential [14]. In particolar, when 
Pi.r (1) = 0, the only solutions are the Bloch states. 
We now consider the case h = h = 0, which gives: 
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In Fig. 2, we present solutions which are symmetric with 
respect to the axis x = 1. These have a linear limit, which 
is simply given by the superposition of two eigenfunctions 
having the same energy but opposite momentum. Obvi- 
ously, in the nonlinear case the superposition principle 
breaks down, so that the existence of such solutions was 
not obvious. 



The density profiles of such solutions are shown in the 
Fig. (3), with different numbers of nodes. 
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FIG. 3. Symmetry broken generalized Bloch states with 
three nodes (a) and with five nodes (b). 

In Fig. 4 we plot the chemical potential of the system 
as a function of nonlincarity. The full lines corresponds 
to zero-current Bloch states, while the dashed and dotted 
lines correspond to generalized Bloch states. Notice that 
when i] ~ —8, the ground state (q = 0, Bloch state) is 
replaced by a symmetry broken Bloch state. 



FIG. 4. Chemical potential as a function of nonlinearity 

Conclusions We have studied the nonlinear Kronig- 
Pcnncy model. This is given by the nonlinear 

Schrodinger equation with a periodic delta-function ex- 
ternal potential. We ve found analytical solutions of zero- 
current states having the same (Bloch states) or a differ- 
ent ( "generalized" Bloch states) periodicity of the poten- 



tial. Nonlinear Bloch states reduce, in the linear limit, 
to the well known eigenfunctions of the linear Kronig- 
Penney model. We have studied the chemical potential 
dependence of such states and compared it with the lin- 
ear Kronig-Pcnney excitation spectrum. 
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